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Abstract 

Structures are often characterized by parameters, such as mass and 
stiffness, that are spatially distributed. Parameter identification of 
distributed structures is subject to many of the difficulties involved 
in the modeling problem, and the choice of the model can greatly affect 
the results of the parameter identification process. Analogously to 
control spillover in the control of distributed-parameter systems, 
identification spillover is shown to exist as well and its effect is to 
degrade the parameter estimates. Moreover, as in modeling by the 
Rayleigh-Ritz method, it is shown that, for a Rayleigh-Ritz type 
identification algorithm, an inclusion principle exists in the 
identification of distributed-parameter systems as well, so that the 
identified natural frequencies approach the actual natural frequencies 
monotonical ly from above. 


Introduction 

Structures are distributed-parameter systems, which implies that 
they are characterized by parameters, such as mass and stiffness, that 
are spatially distributed. The motion of distributed-parameter 
structures is governed by partial differental equations (pde s), in 
which the parameters appear in the form of space-dependent coefficient 
functions. For the most part, these pde's do not admit closed-form 
solutions, so that one must be content with approximate solutions. To 
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obtain an approximate solution, it is necessary to resort to spatial 
discretization, which amounts to representing the motion of the 
structure by a linear combination of admissible functions depending on 
the spatial variables multiplied by time-dependent generalized 
coordinates (Ref. 1). Certain integrations over the structure eliminate 
the dependence on the spatial variables. The net result is the 
transformation of the pde's into a set of ordinary differential 
equations (ode's), in which the parameters enter into the coefficient 
matrices. The process of obtaining the set of ode's is equivalent to 
deriving a discretized model designed to approximate the distributed 
model of the actual structure. How well the discrete model is capable 
of representing the distributed model depends on the number and type of 
admissible functions used in the discretization process. Indeed, the 
rate of convergence to the exact solution depends on the number of 
admissible functions and, perhaps to a larger extent, on the type of 
admissible functions used (Ref. 2). The discretized model can be used 
to compute a finite number of lower modes of vibration. It is 
symptomatic of such approximations that an even smaller number of 
computed lower natural frequencies than the number of degrees of freedom 
of the discrete model are accurate representations of the actual natural 
frequencies (Ref. 1). The above discussion pertains to structural 
modeling on the basis of given distributed parameters. The discrete 
model can be used not only to compute the modes of vibration but also to 
derive the system response. The problem of deriving the response to 
given excitations when the system characteristics are known, albeit only 
approximately, can be regarded as a direct problem. 
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On occasions, the system parameters are not known and the object is 
to infer them from measured system response to known excitations. This 
represents an inverse problem, more commonly known as parameter identi- 
fication, or parameter estimation. Clearly, in the case of distributed 
structures, the object is to identify parameter distributions. The 
implication is that the parameters are distributed nonuniformly, because 
otherwise the problem would be almost trivial, reducing to the 
identification of a single number per parameter distribution. But, 
nonuniform parameter distributions are precisely those preventing 
closed-form solutions and demanding discrete models. Hence, if the 
object is to identify nonuniform parameter distributions, then one must 
expect at the very least the same kind of idiosyncrasies encountered in 
mere modeling, so that one must proceed with extreme caution. 

Commonly used discretization procedures, for modeling as well as 
for parameter identification, are the classical Rayleigh-Rltz method, 
the finite element method, the Galerkin method, etc. (Ref. 1). They all 
end up characterizing the system parameters by means of matrices, such 
as the mass and stiffness matrices. 8ut, this is the very same process 
that proved incapable of giving wholly accurate eigensolutions. Indeed, 
it is typical of the above techniques that less than one half of the 
computed natural frequencies, and even fewer natural modes, are 
sufficiently accurate (Ref. 1). Hence, it is unreasonable to expect 
that a process identifying parameter matrices can yield more accurate 
results than those obtained by a corresponding process in mere modeling. 

There are several factors that call for a cautious approach to 
parameter identification in distributed structures. In particular, one 
must interpret and use the results of the identification process 
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judiciously. The mass and stiffness matrices are not unique for a 
structure, and in fact they do not even represent physical quantities. 
Indeed, the dimension of the matrices depends on the number of admis- 
sible functions used and the entries of the matrices depend on the 
nature of the admissible functions. These two factors are in fact 
related, as a wise choice of admissible functions can give superior 
results with a smaller number of admissible functions. In particular, 
one must make sure that the admissible functions are such that all the 
boundary conditions can be satisfied by finite linear combinations of 
these functions (Ref. 2). The mass and stiffness matrices themselves 
have no physical meaning. For distributed structures, the quantities 
having physical meaning are the mass and stiffness distributions. 

Perhaps the idea can be best illustrated by using the analogy with the 
construction of a building, in which the mass and stiffness matrices 
represent the scaffolding and the mass and stiffness distributions 
represent the building itself. Hence, an identification process must 
not have as its goal the identification of the mass and stiffness 
matrices, but the identification of the mass and stiffness distributions 
(Ref. 3). At times, the identification of natural frequencies and modes 
can serve as an intermediate step (Refs. 4 and 5). Another factor 
affecting the quality of the results is the design of the experiment. 

In particular, one must ensure that the actuators and sensors are at 
significant points of the structure and that they are in sufficient 
number. Of course, noisy actuators and sensors can cause difficulties 
and even lead to instability in the identification process (Ref. 6). 

In modeling, the question arises as to the effect of truncation on 
the system eigenvalues. Specifically, what happens to the system 
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eigenvalues when the number of admissible functions is increased (or 
decreased) in the modeling process? The answer to this question lies in 
the inclusion principle (Ref. 1), which can be used to verify the 
accuracy of a model. Because parameter identification in distributed 
structures depends heavily on modeling, the question can be asked 
whether an inclusion principle exists in parameter identification. 

Practical limitations dictate that control of distributed 
structures be designed and implemented on the basis of truncated models, 
whereby only a subset of the modes of the structure are actively 
controlled. The excitation of the residual (uncontrolled) modes caused 
by the finite-dimensional controller is known as control spillover and 
can degrade the system performance. The contamination of the sensor 
measurements by the residual modes is known as observation spillover, 
and the combination of control and observation spillover can destablize 
the system (Ref. 7). An analogous effect exists in the identification 
problem. Because distributed structures must be modeled by discrete 
systems, the contamination from the residual (unmodeled) modes leads to 
identification spillover. Identification spillover can cause errors in 
the parameter estimates. 

This paper stresses the importance of judicious modeling for 
parameter identification. Distributed-parameter systems can be 
approximated by a variety of models. The accuracy of these models 
depends on the type and number of admissible functions used in the 
discretization process. As shown in the numerical example, the accuracy 
of the identified parameters depends on the type and number of 
admissible functions. In parameter identification, the system 
parameters can be estimated directly using the system response to known 
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input (Refs. 3 and 6). In this paper, the modal quantities are 
identified first, and then the system parameters are identified on the 
basis of the estimated modal quantities (Ref. 4). 

Modal identification of distributed structures can be performed 
using a variational approach (Ref. 5). A ratio of the maximum potential 
energy to the reference kinetic energy, known as the Rayleigh quotient, 
can be used to generate the eigenvalues of the model. Based on the 
inclusion principle, we conclude that, when the order of the model 
increases, the eigenvalues of the model approach the eigenvalues of the 
distributed structure monotonically from above. The variational modal 
identification technique employs the temporal and spatial properties of 
distributed structures to generate a pseudo-Rayleigh quotient, in which 
stationary values of the quotient approximate the eigenvalues of the 
actual system. Moreover, as the order of the identified model is 
increased, the identified eigenvalues approach the actual eigenvalues 
monotonically from above. A numerical example demonstrates that the 
inclusion principle holds for the variational modal identification 
technique. 

Equation of Motion 

We consider the case in which the linear motion of the distributed 
structure is governed by the partial differential equation 

£u(P,t) + m(P)u(P,t) = f(P,t) , P € 0 (1) 

where £ is a self-adjoint positive definite differential operator of 
order 2p, m is the mass density, u is the displacement at the spatial 
position P at time t, f is the external force density and D is the 
domain of extension of the system. The displacement u satisfies the 
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boundary conditions 


B^u(P,t) - 0, i - 1,2,. . • ,P » P S (2) 

where B. are homogeneous differential operators of maximum order 2p-l 
and S is the boundary of 0. The boundary conditions are either geo- 
metric, in which case the order of B. is smaller than p, or natural, in 
which case the order of B. is smaller than 2p. Functions that are p 
times differentiable and satisfy the geometric boundary conditions are 
called admissible functions. Functions that are 2p times differentiable 
and satisfy all the boundary conditions are called comparison functions. 
Associated with Eq. (1), we have the eigenvalue problem 

MP) = xm(PMP) ( 3 ) 

where <s(P) are functions satisfying the boundary conditions, Eq. (2). 

The solution to Eq. (3) consists of a denumerably infinite set of real 
eigenfunctions a> r (P) and associated real positive eigenvalues 

(r = 1,2,...). For convenience, we order the eigenvalues so that 
X 1 ~ v 2 “ **' * The eigenfunctions are orthogonal and they can be 

normalized so as to satisfy the orthonormality conditions 

j* m(PU (P)*_(P)dD = S , J * r (py:* s (P)dD = x r 5 r$ (4) 

D 0 

where & is the Kronecker delta function. From the expansion theorem 
r s 

(Ref. 1), we can express the displacement u as a linear combination of 
the eigenfunctions, or 

u(P,t) = l ♦ (P)u r (t) ( 5 ) 

r=l 

where u r .(t) are modal coordinates. Introducing Eq. (5) into Eq. (1), 
considering Eqs. (4) and following the usual steps, we obtain the 
infinite set of second-order ordinary differential modal equations 
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( 6 ) 


u r (t) + uyJ r (t) = f r (t), r - 1,2,... 

2 

where x r = u£, in which uy represent the natural frequencies of 
oscillation, and 

f r (t) = / V (P)f(P,t)dD, r = 1,2, .. . (7) 

0 

are modal forces. 

In the case of free vibration, Eqs. (6) reduce to 

u r (t) + uyu r (t) =0, r = 1 , 2 ,... (8) 

Equations (8) are independent, so that distinct modes of vibration are 
uncorrelated, or 

, T 

< u r (t), u $ (t) > = lim f f u r (t)u (t) dt = Q r « , r.s = 1,2,... 

T “ ° (9) 

where the symbol < , > represents a temporal inner product and Q r (r = 
1,2,...) are positive constants. Using the solution of Eqs. (8), it can 
be shown that 

< u r (t), u $ (t) > = uj;q r 5 rs , r,s = 1,2,... (10) 

The modal coordinates u r (t) can be obtained from the system response by 
using the modal filters (Ref. 8) 

u (t) = f m(P)<t> (P)u(P,t)dD , r = 1,2,... (11) 

D 

The eigenvalue problem given by Eqs. (2) and (3) can be replaced by 
a variational problem consisting of determining the stationary values of 
the Rayleigh quotient 

R(m(P)) = — »( p )l (12) 

f m(P)/(P)dD 
0 

where ^(P) is a trial function and [ , ] is an energy inner product 
(Ref. 1), the latter being obtained through integrations by parts 
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of J *(PVU(P)dD. Note that the energy inner product I , 1 Is 
symmetric because the stiffness operator 4. in Eq. (1) is seif- 
adjoint. Moreover, the energy inner product [u, u| is proportional to 
the potential energy of the system (Ref. 1). Expressing the trial 
function in Eq. (12) as a linear combination of the eigenfunctions, or 

*(P) = l c r » r (P) (13) 

r= I 

and considering the spatial orthonormality conditions, Eqs. (4), Eq. 


(12) becomes 

R ( c ^ , Cg » • * • ) 



2 2 
Vr 



(14) 


We note that the stationary values of the Rayleigh quotient R are 


identical to the system eigenvalues and occur every time the trial 


function *(P) is identical to an eigenfunction (Ref. 1). 

Closed-form solutions to Eq. (1) exist in a few simple cases only, 
namely cases in which the parameter distributions are uniform. For more 
complex structures, an approximate solution must be sought. To this 
end, the partial differential equation must be discretized in space, so 
that an approximate solution is obtained by means of a finite set of 
ordinary differential equations. Commonly used discretization 
procedures are the classical Rayleigh-Ritz method, the finite element 
method and the Galerkin method. The difference between the Rayleigh- 
Ritz and Galerkin method lies in the choice of trial functions, in the 
sense that Galerkin trial functions must be comparison functions, 
satisfying all the boundary conditions. If the Rayleigh-Ritz method 
employs the variational formulation given by Eq. (12), then the trial 
functions can be merely admissible functions, satisfying the geometric 
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boundary conditions alone. Of course, the Rayleigh-Ritz method is 
applicable only to self-adjoint systems, whereas Galerkin's method can 
be used for non-self-adjoint systems as well (Ref. 1). The finite 
element method is a variant of the classical Rayleigh-Ritz method, and 
the trial functions are local admissible functions, as opposed to the 
trial functions in the classical Rayleigh-Ritz method, which are global 
admissible functions. 

Approximate Methods. The Ra,yleiqh-Ri tz Method 

Consider the admissible functions i> r (P) (r = 1*2 n) satisfying 

the geometric boundary conditions and express the displacement u(P,t) as 
a linear combination of the admissible functions 4y(P) multiplied by 
time-dependent generalized coordinates d r (t) as follows: 

u (n) (P,t) = l * r (P)q r (t) - S T t p >9 (t > (15) 

r= 1 

where u^ n ^(P,t) is the n th -order approximation of u(P,t), 

'O(P) = Ui(P) 0> 2 < p ) ••• ^(P) 1 T is an n-vector of admissible functions 
and q(t) = (q^t) q 2 (t) ... q p (t) 1 T is an n-vector of generalized 
coordinates. As n approaches infinity, the Rayleigh-Ritz solut- 
ion i/ n) (P,t) converges to the exact solution u(P,t), provided the trial 
functions * (P) (r = 1,2,...) are complete in energy (Ref. 1). 
Substituting Eq. (15) into Eq. (1), multiplying on the left by *(P) and 
integrating over the domain 0 of the structure, we obtain the n-degree- 
of-freedom discretized model 

. K (n) g(t) » F< n) (t) ( l6 > 

where 

M (n) , j m(P)*(P)* T (P) dO, = U(P), h> T (P)1 (17a, b) 

D 

are n x n mass and stiffness matrices, respectively, and 
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(18) 


F (n) (t) = j" #(P)f(P.t) dO 
0 " 

is an n-dimensional generalized force vector. 

An attractive feature of the Rayleigh-Ritz method is that, by 
increasing the number of degrees of freedom n, the previously calculated 
mass and stiffness coefficients do not change, so that one need only 
calculate an additional row (or column) to obtain the updated mass and 
stiffness matrices. For an (n + 1) -degree-of-freedom discretized model, 
the mass and stiffness matrices have the form 


' m<"> ; 

K (n+1) . 

k<"> i 

x x x_ 


XXX 


The eigenvalue problems associated with the discrete 


(19a, b) 

models are given by 


(n) M (n) (n) = *(")>), r = 1,2,. ...n 
r -r 

v (n+l) M (n+l) x (n+l) s K (n+1) x (n+l) f r = 1,2 n + 1 


(20a) 
(20b) 

~r -r 

The orthonormality conditions require that the solutions to Eq. (20a) be 
normalized so that 

? (n) T H (n) ? (n) . ^ ;(">>>;<"> - »< n) « rs ( 21a - b > 

From the inclusion principle (Ref. 1), the relationship between the 

eigenvalues k< n) (r • 1,2 n) of the n-degree-of-freedom model and 

the eigenvalues x< n+1) (r • 1,2 n * 1) of the (n * l)-degree-of- 

freedom model is given by 


(-1) S xS"> S 4 ntl) s 4 n) S ... * 4"’ 4 ‘n^’ 


( 22 ) 


i. - "'l ■ A 2 - a 2 “ ~ 'n n+1 

As the number of degrees of freedom n in the discretized model is 
increased, the approximate eigenvalues decrease monotonical ly and 
approach the actual eigenvalues of the system asymptotically from above 

(Ref. 1). 
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For complex structures, the system parameters, which in the case at 
hand are the mass and stiffness distributions, are complicated functions 
of the spatial variables. Because closed-form solutions to Eq. (1) 
exist in simple cases only, almost always involving uniform parameter 
distributions, we must obtain an approximate solution using Eq. (16). 
This implies that the use of Eq. (1) for identification purposes is not 
practical, so that an exact parameter identification is not feasible. 
The inverse problem of identifying the unknown parameters represents a 
compounded version of the direct or modeling problem, where in the 
latter the parameters are known. We must expect at least the same 
idiosyncrasies in parameter identification as those encountered in 
modeling. A suitable choice of the number and type of trial functions 
used in modeling to obtain Eq. (16) must be determined prior to any 
parameter identification process. 


A Rayleigh-Ritz Type Paramet er Expansion 
As in the Rayleigh-Ritz method, we expand the parameter 
distributions in terms of a set of known functions multi pled by 
undetermined coefficients. It is assumed that accurate representations 


of the mass and stiffness distributions are given by (Ref. 4) 


m(P) - l W(P), kC) * l. 3 r k r (P) 
r=l r= l 


(23a, b) 


where m (P) and k f (P) are functions from complete sets and o r and 8 r 
are undetermined coefficients playing the role of unknown parameters in 
the identification process. The functions m r (P) and k r (P) can be global 
or local functions. Using Eqs. (17) and (23), it can be shown that the 
coefficient matrices in Eq. (16) can be written in the form (Ref. 3) 


» = 


9 


(n) „(n) 


■ l VV • 10 

r=l 


I 3 r r 
r=l 


(n) 


(24a, b) 
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where 


= J m (P)-4»(P)'4» T (P) dD . K r = M p ). /( p )l r (25a, b) 

r D r ~ 

Equations (24) are not new. In fact. Ref. 9 uses submatrix scaling to 
define distributed parameters where the mass and stiffness matrices are 
expanded into a linear series of submatrices as given by Eq. (24). One 
must proceed with caution, however, as Eqs. (15) and (23) represent 
approximations to the displacement profile and distributed parameters, 
respectively. Indeed, the discretized model given by Eq. (16) was 
obtained as an approximate solution to Eq. (1). Moreover, the dis- 
cretized representation given by Eq. (16) is not unique as the entries 
in the mass and stiffness matrices depend on the type of admissible 
functions used in Eqs. (15) and (25). The accuracy of the discretized 
model depends on the type and number of admissible functions used (Ref. 
2). It follows that the accuracy of the estimated parameters in the 
identification process also depends on the type and number of admissible 
functions used. 


Parameter Identification 

In this section, we present an algorithm for identifying parameter 
distributions in distributed structures. The algorithm represents an 
iterative procedure in which the parameter distributions are updated 
using a set of natural frequencies identified from the measured response 
of the distributed structure. In the next section, we present a method 
for identifying the natural frequencies from the system response. In 
the ensuing discussion, we refer to the identified natural frequencies 
of the distributed structure as the measured natural frequencies. We 
adopt an incremental approach in which a vector u = u > ^ ••• 1 * 
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identified natural frequencies of the actual distributed structure, 
playing the role of a measurement vector, is obtained and then used to 
update a parameter vector given by p = °2 • •• <*g a 2 • •• s h l 

The iteration process is based on the incremental relations 


^ = » - w (n) , Aw = (^Iap , AP = P - P 0 (26a, b,c) 

j 

To begin the parameter identification, we postulate a parameter 
vector p Q , compute the corresponding mass and stiffness matrices using 
Eqs. (24) and (25) and solve for the associated natural frequencies 
u (n) = [/TnT /T"! _ XjnT] US ing Eq. (20a). Then, we use the 

measurement vector w to compute Aw by using Eq. (26a). The Jacobian 
matrix in Eq. (26b) can be computed by using the orthogonality 
conditions, Eqs. (21). The entries are 


TpT 


i 

2w 


,(n) T 3K^_ 

; 1 ap« 


,(n) . 


a 2 x ( n > 

u i 


aM 

ap. 


(n) 


- X 


(n) 


i = 1,2,... ,f j = 1,2,. ...9 


(27) 


where the eigenvectors x( n ^ are computed using Eq. (20a) and the entries 
of the mass and stiffness sensitivity matrices are given by 


aM 


(n) 


ap 

aM 


j 

( n ) 


ap 


J 


= 0, 


aK 


(n) 


ap 


- 0, j = 


aK (n) _ „(n) 


(28) 


ap. 


= k . , j = g + i.g + 2, •••.9 + h 


'j -'J 

Then, the increment Ag can be computed by means of Eq. (26b). For the 
case in which the number f of measured natural frequencies is greater 
than the number of parameters g + h, a least-squares solution of Eq. 
(26b) yields 
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ap = 


3w ^ 


3 cjj ^ 

W. 


i-l 


3oj • 

3fT 


Ad) 


(29) 


Finally, from Eq. (26c), we obtained the improved parameter vector 


p = P Q + Ap 


(30) 


The procedure is repeated using Eqs. (26a), (27), (29) and (30) until 
convergence is achieved. 


An Inclusion Principle in Parameter Iden tification 
In the preceding section, the estimated parameters were updated 
using the identified natural frequencies of the distributed structure. 
The objective is to identify the natural frequencies u> r (r = l,2,...,f) 
from the free response. To this end, we define a pseudo-Rayleigh 
quotient suitable for modal identification in the form (Ref. 5) 


RU(P)) 


< w(t) , w(t) >_ 

< w(t), w(t) > 


(3D 


where 


w(t) = I c(P)u(P.t)dD 


(32) 


in which C (P) is an admissible function, i.e., a function satisfying the 
geometric boundary conditions. We can express the admissible 
function c(P) as 


c(P) - I a r m(P)^> r (P) 

r=l 


(33) 


where n r (P) are eigenfunctions. Introducing Eqs. (5) and (33) into Eq. 
(32) and the result into Eq. (31), considering the orthonormality 
conditions, Eqs. (4), and the temporal correlation properties, Eqs. (9) 

and (10), we obtain 


R ( * • • • ) 


? a 2 Q J 

L. a r w r r 

r= 1 

* o 

l. 

r=l 


(34) 
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Note that Eq. (34) Is identical to Eq. (14) with c f • a r /Q r , so that 
stationary values of the pseudo-Rayleigh quotient, Eq. (31), are 
Identical to the eigenvalues = w P and they occur when c(P) is equal 
to an eigenfunction multiplied by the mass distribution. 

In practice, it is only possible to identify a finite number of 
natural frequencies and associated modes of vibration. To this end, we 
express the trial function c(P) as a finite linear combination of the 
admissible functions i» r (P) as follows: 

t ( ">(P) ■ j 1 » r * r< P > <35 ’ 

where ^ m V) is the m th -order approximation of c(P) and where v r are 
undetermined coefficients. Introducing Eqs. (15) and (35) into Eq. 


(32), we have 

w(t) = l ? v^q^t), Vi ■ *r( p >*1 (p)<i0 

r=l i=l u 

Introducing Eqs. (36) into Eq. (31), we obtain 

y y „<»> v v 

, r.l s- 1 fS M 

P( v ^» v 2» * • * ,v n^ m m 


(36a, b) 


(37) 


where 


I, l "rs’ Vs 

r=l s=l 


n n 


m<™> - .<;> - I l *r,*,j« < ’i< t >* < ’j (t) r 

3 i=l j=l 


,s, - 1 , 2 , ... ,m 

(38) 


k (m) = k (m) = l J Vi 4 'sj < >J r,S ’ = 

Determining the statiJwy values of Rayleigh's quotient, Eq. (37), can 

be replaced by the eigenvalue problem (Ref. 1) 

x (m) M (m) v (m) = K (m)^(m) (39) 

where the solution to Eq. (39) represents the m th -order approximation to 
the stationary values of the pseudo-Rayleigh quotient, Eq. (31), and 
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where M^ m ), K( m ) and v^ have corresponding entries [■JlJ'I. (kj^l 
and v< m >, respectively. The m th -order approximate eigenfunctions can be 
computed from the eigenvectors v^ using 

*M(P) - - T (P)^ m) ( 40 > 

The solution to Eqs. (39) and (40) represent the identified mth-order 
approximation to the eigenvalues and eigenfunctions of the actual 
system, in the same way as the solution to Eq. (20a) represents the nth- 
order approximation to Eq. (3). 

When m = n in Eqs. (38), the orders of the model and identified 
model are equivalent, so that in this case identification spillover does 
not occur. For the model to be an accurate representation of a 
distributed structure, the number of degrees of freedom n can be quite 
large, so that in practice the number m of identified natural 
frequencies and mode shapes can be much smaller than n. In the case in 
which m < n, identification spillover is present. Note that, by 
increasing m in Eqs. (38), previously calculated entries in 
and do not change. Hence, as in modeling (Ref. 1), the identified 
eigenvalues satisfy the inclusion principle as given by inequalities 
( 22 ). 

Numerical Example 

As an illustration, we consider a rod in axial vibration with a 
spring of spring constant k attached to the free end as shown in Fig. 

1. For this case, the stiffness operator L in Eq. (1) has order 2 (p = 
1) and is given by 

[woo y (41) 

where EA(x) is the axial stiffness. The boundary conditions are given 
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by 


8^=1 at x = 0, B 1 = EA + k at x L 


Moreover, the energy inner product in Eq. (12) has the form 


(u, u] = f EA 
0 


1 ) 


’dx + ku^(L.t) 


(42) 


(43) 


The mass and stiffness distributions for the actual model are given by 

2 


/ \ 6 m r, ifX 

(x) = tt l 1 - ?ir 


l^Xi2i 


2^L 


EA(x) 


6EA r i 1 < x 1 


U - jlr) I 


(44a, b) 


As an illustration, we consider the case where the distributions are 
known with unknown scaling such that g = h = 1 and that 

yx) = | [i - \'i )], m*) - f ii - ) ( r )' in Eqs - (23) - Hence> the 

problem reduces to identifying three parameters, namely, m, EA and k. 
Moreover, because free vibration data is used, to identify the 
parameters uniquely, one of the parameters must be known. For 
convenience, we chose m = 1 kg/m. 

The first objective is to show that accurate modeling has a direct 
influence on the accuracy of the identification results. The system 
does not have a closed-form solution so that we must resort to spatial 
discretization. To this end, we consider two sets of trial functions 
* (P) (r « 1,2 n) to be used in Eq. (15). The first set consists of 

the admissible functions 

, . . (2r - l)wx „ _ i ? n (45) 

4y(x) * sin -» — 21 — — » r " 

We note that these functions satisfy the geometric boundary condition u 
= 0 at x = 0. They represent the eigenfunctions to a related problem, 
i.e., they are the eigenfunctions of a uniform rod with no spring at the 
free end. The second set of trial functions consists of the functions 

(46) 


5 r (x) = sin 


rnx 

2L ’ 


r = 1,2, ... ,n 
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These functions are also admissible functions and they represent quasi- 
comparison functions for this example (Ref. 2). Note that a finite 
linear combination of the quasi -comparison functions can satisfy the 
natural boundary condition (i.e., = EA + k at x = L) , whereas the 

first set of admissible functions, although complete in energy, can only 
satisfy the natural boundary condition when n approaches infinity. 

To refine the system parameters EA and k, we used the first three 
natural frequencies from the actual model so that f = 3 in Eq. (28). As 
an initially postulated model, we used the parameters EA = L.5 N and k = 
1.5 N/m, while the actual model parameters are EA = 1 N and k = 1 N/m. 

To check the algorithm, several other cases involving different 
initially postulated parameters were tested. The results are not 
presented here, but for reasonable starting values the identification 
results were insensitive to the initially postulated parameters and the 
algorithm converged within five iterations. In addition to using 
ordinary admissible functions and quasi -comparison functions for 
modeling, different numbers of degrees of freedom were used. Table 1 
displays the results. As expected, the identified parameters improve as 
the number n of degrees of freedom increases. Indeed, as the number n 
of admissible functions increases, the model becomes more representative 
of the actual system, resulting in improved identified parameters. The 
model with the admissible functions can predict the value for EA very 
well, but the identification of the parameter k is rather poor because a 
finite number of admissible functions cannot satisfy the natural 
boundary condition. On the other hand, the parameter identification 
based on the model using quasi -comparison functions is extremely good 
for n > 4, as both parameters are identified exactly. 
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The second objective is to illustrate the inclusion principle in 
parameter identification and demonstrate the effects of identification 
spillover. To this end, we consider the rod in axial vibration 
described previously. To identify the natural frequencies from the free 
response, it was necessary to simulate the response of the actual 
system. The model for this purpose was obtained by the Rayleigh-Ritz 
method using n = 8 admissible functions in Eq. (15). The admissible 
functions used in the discretization were the quasi -compari son functions 
described by Eq. (46). Different order models were identified using the 
free response, where m = 4,5,6 and 7 degrees of freedom were considered 
in Eqs. (39). With the quasi-comparison functions as admissible 
functions, Eq. (36b) yields 



We note that, to compute q^(t) and q - ( t ) (i = 1,2 8) in Eq. (38), 

the model was excited using an initial impulse of magnitude 1 N-s acting 
at x = L. To identify the natural frequencies, the solution of the 
eigenvalue problem given by Eq. (39) was obtained and the results are 
presented in Table 2. Note that identification spillover has the effect 
of degrading the identified results. Indeed, the spillover from the 
unmodeled degrees of freedom has the effect of increasing the identified 
natural frequencies. Moreover , as the order of the identified model 
increases, the estimated natural frequencies decrease monotonically and 
approach the actual natural frequencies asymptotically from above. It 
is also obvious from Table 2 that the computed eigenvalues satisfy the 
inclusion principle, as stated by inequalities (22). 
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Conclusions 


Parameter identification in distributed structures is ordinarily 
based on discretized models defined by mass and stiffness matrices. 

Mass and stiffness matrices, however, have no physical meaninq for 
distributed structures, whereas mass and stiffness distributions do. 
Indeed, the mass and stiffness distributions are the quantities 
describing the physical characteristics of the structure. Moreover, 
mass and stiffness matrices are not unique for a given structure and 
their entries and dimensions depend on the number and type of trial 
functions used in the modeling process. Hence, the object of a 
parameter identification technique should be to identify physical 
properties and not mass and stiffness matrices. 

Parameter identification in structures represents a compounded 
version of the modeling problem in structures. Indeed, as the numerical 
example indicates, the choice of the model can greatly affect the 
results of the parameter identification process. The choice of the type 
and number of admissible functions is at least as important in parameter 
identification as it is in the problem of modeling distributed 
structures. A wise choice of the type and number of the admissible 
functions, can improve the results of the parameter identification 

process. 

Practical limitations dictate that the control of distributed- 
parameter systems be designed on the basis of finite-order models. The 
excitation of the unmodeled modes by finite-dimensional controller can 
degrade the system performance. As in the control problem, 
identification spillover exists and it has the effect of degrading the 
estimates of the parameters. It is shown that an inclusion principle 
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exists in the identification of distributed-parameter structures, so 
that the identified natural frequencies approach the actual natural 
frequencies monotonically from above. 
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Table 1. Identified Parameters Using Admissible Functions (AF) and 
Quasi -Comparison Functions (QCF) for Oifferent Oegrees of 


Freedom. 

AF 


n 

EA 

k 

% RMS 
ERROR 

4 

1.0016 

0.91705 

5.87 

5 

1.0010 

0.93335 

4.71 

6 

1.0008 

0.94474 

3.91 

7 

1.0006 

0.95286 

3.33 

8 

1.0005 

0.95208 

2.89 


QCF 


EA 

k 

% RMS 
ERROR 

0.98809 

1.0564 

4.08 

0.99996 

0.99977 

0.02 

0.99996 

0.99984 

0.01 

0.99999 

0.99987 

0.01 

1.00000 

1.00000 

0.00 


Table 2 Identified Natural Frequencies for Different Oegrees of Freedom. 

Identified Natural Natural Frequencies of 

Frequencies (rad/s) Discrete Model (rad/s) 


m = 4 

m = 5 

m = 6 

m = 7 


2.388 

2.250 

2.249 

2.228 

2.210 

15.083 

5.225 

5.223 

5.133 

5.083 

21.886 

17.165 

8.822 

8.173 

8.087 

57.160 

23.860 

17.198 

12.860 

11.148 


63.161 

23.918 

17.308 

14.237 



63.351 

24.155 

17.535 




64.156 

24.388 





64.946 
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